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A simple picture for the spectrum of the one-dimensional Hubbard model is presented using 
a classification of the eigenstates based on an intuitive bound-state Bethe-Ansatz approach. This 
approach allows us to prove a "string hypothesis" for complex momenta and derive an exact for- 
mulation of the Bethe-Ansatz equations including all states. Among other things we show that 
all gapped eigenstates have the Bethe Ansatz form, contrary to assertions in the literatureu. The 
simplest excitations in the upper Hubbard band are computed: we find an unusual dispersion close 
to half-filling. 
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I. INTRODUCTION 

In a recent articlei we have introduced a new approach to the derivation of the Bethe Ansatz equations of models 
with bound states. By bound states we mean that the wave functions contain complex momenta kj which need to be 
properly handled in the infinite volume limit. Our approach avoids postulating a "string hypothesis" for the momenta 
and allows a systematic and simple procedure to derive all eigenstates. 

In our approach the Bethe Ansatz is formulated for bound states (and bound states of bound states). These 
composites are formed on the infinite line and are incorporated systematically as building blocks of the wave functions 
defined on a finite ring. The composites correspond to poles of the respective S-matrices thus guaranteeing their 
stability under scattering with other particles. This necessitates the introduction of appropriate boundary conditions 
- Composite Boundary Conditions (CBC) - which respect the construction and allow the formation of exact strings. 
We conjecture that the CBC allow a symmetry - presumably Yangian - to be manifest already at finite volume, a 
symmetry that is violated by the conventional Periodic Boundary Conditions, hence appearing in the latter case only 
in the infinite volume limit. We shall refer to this approach as the Bethe Ansatz for Composites (BAC). 

The approach is general and applies to any model where complex momenta appear: the Hubbard model, the t- 
J model, the Anderson model and the multichannel Kondo model among many others. We shall discuss here the 
Hubbard model in detail and leave the treatment of the other models for later. 

The article is organized as follows: 

In section II we present the model and construct the class of bound state solutions to the Schrodinger equation of 
the model and show how to incorporate them as new ansatz functions besides the usual plane waves into the BA. We 
introduce appropriate boundary conditions (the CBC) to define the model on a finite configuration space, and deduce 
the BAC equations. A particular class of solutions of the BAC equations is found to underlie the-ry-pairing and the 
charge SU(2) symmetry group. Next, we clarify the connection with Takahashi's string hypothesise!, which, we argue, 
leads to inconsistencies when finite volume corrections are taken into account. 

In section III we study the bound state excitation spectrum both in the repulsive and attractive models. In the 
repulsive case, holding the number of electrons fixed, we find a three parameter excitation, consisting of two gapless 
holons (spinless, carrying charge -e) and a bound state residing in the upper Hubbard band (spinless, carrying charge 
2e). The latter is in fact an anti-bound state, corresponding to the formation of a bound pair of electrons with positive 
binding energy dressed by its interaction with the sea electrons. In the attractive case a dual picture emerges: it 
corresponds to breaking one of the pairs forming the ground state. The resulting excitation consists of two gapless 
(renormalized) electrons each carrying charge e and spin 1/2 and a gapped spinless excitation with charge -2e. The 
gapped excitation is an independent mode only away from half-filling and we give numerical results for its dispersion. 

In section IV we summarize our results. Appendix A contains a detailed demonstration of the stability of the 
two-particle bound state in the presence of other particles with composite boundary conditions. 



1 



II. DERIVATION OF THE BAC EQUATIONS 



The one-dimensional Hubbard hamiltonian is given by 

oo 

h= ^2 + h - c -) + Un r,i n i,i- (!) 

i— — oo 

The hamiltonian was diagonalized by Lieb and WuO in the sector with ./V particles on a finite ring of length L, with 
periodic boundary conditions (PBC) imposed on the wavefunction of the iV electrons, 

F(x-L,...,Xi,...,x N )=F(xi,...,Xi + L,...,x N ) V i. (2) 

The resulting eigenfunction isjmrameterized, for total spin S — h(N— 2M), by N momenta kj and M spin rapidities 
A 7 , satisfying the BA equations™], 

Ag hill fi,j -f f-j 



n 



A 7 - A5 + it 



n 



sin 



A-v — sin kj + i 1 



(4) 



where u — U/t. 

The energy and momentum of the state are then given by 

E = -It ^2 cos k j p = X! fc J ■ 

The eigenstates of the hamiltonian correspond to the various solutions of eqns (^,^) . Real as well as complex solutions 
need to be considered to obtain a complete spectrum of states. 

We shall argue, however, that complex spin rapidities {A 7 } and complex charge momenta {kj} have a different 
character although both of them were treated in much the same way in the literature. The former type is associated 
with the spin degrees of freedom and describes kink/anti-kink bound states. As such they are a many-body phe- 
nomenon and make their appearance through complex conjugate solutions of the BA equations for spin rapidities. 
The existence of these solutions is postulated in the A string-hypothesis about the form of ths,BA solutions for large 
number of particles N. Although very plausible, this hypothesis remains unproven up to nowO. 

The latter type - on which we concentrate in this article - is associated with the charge degrees of freedom. In 
analogy to the spin sector, they were assumed to correspond to complex conjugate solutions of the Lieb - Wu equations 
for charge rapidities. We will show that they have to be treated differently because they are not due to a many-body 
effect but can be identified with the elementary bound states of the hamiltonian, present already in the few-particle 
sector. They have to be incorporated into the Bethe ansatz ab initio and lead to a set of BA equations, which we call 
the Bethe ansatz for composites (BAC). 

The charge complex momenta were also conjectured to form strings, the so called k — A strings, and Bethe- Ansatz 
equations were derived based on this hypothesisu. In the thermodynamic limit the BAC equations turn out to be 
equivalent with those derived using the string-hypothesis for the charge rapidities. This finding does not constitute, 
however, a proof of the string-hypothesis for a finite system. On the contrary, we shall argue in section II D. that 
the string-hypothesis for charge bound states leads to an over- constrained set of equations, which has in general no 
solutions for sufficiently large but finite system size L. 



We begin by discussing the two string hypotheses separately. To introduce "A - strings" one assumes complex- 
conjugate pairs of rapidities = Ao ± i\, and introducing them into eq (^) one concludes, 

X = \ + 0(e-* N ) (5) 

with k > 0. More generally, m spin rapidities Aj are grouped together to form a A-string of length m: 

u 

Aj = A + i(m + 1 - 2j)— j = l... m. (6) 
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The solutions with complex A to the BA equations are driven to the string position (g) in the limit N — > oo, 
corresponding to a many particle effect. 

The standard classification emerges in terms of spinous - chargeless spin-1/2 - objects. The ground state is a 
singlet, 5 = 0, described by a solution with M = N/2 real A rapidities. The simplest spin excitation corresponds to 
M = N/2 — 1 rapidities, hence S = 1. It is found to be a two parameter state and can be interpreted naturally as 
a symmetric combination of two elementary excitations, each carrying spin 1/2, the spinon states. To confirm this 
picture one needs to show that there is also a state with the spinons combined antisymmetrically to form an excited 
singlet. This can be done by introducing a A-string of length 2 so that again M = N/2 corresponding to a singlet, 
S = Op. This structure of excitations is universal to SU{2) invariant models. It was first observed in the GrossJiJeveu 
model], and subsequently rediscovered in several other models: in the Kondo modelEI, the Heisenberg modelE], the 
Hubbard modeltil or the t-J modeli3. Also the spinon - spinon S-matrix, when expressed in terms of rapidities, is 
identical in these models, and was originally calculated from the Bethe Ansatz in [ ||. The fact that quantities such 
as the quantum numbers of the spinons and their scattering phase shifts are the same in the different models is due 
to the circumstance that in all of them the interaction has the form of a spin exchange. Thus the spin sector of the 
various models obeys the same BA equations up to shifts specific to each. As phase shifts and quantum numbers are 
deduced from simple counting arguments that are independent of those shifts they are identical in spite of a wide 
variation in dynamics, see [ ||] for a detailed discussion. 

In addition-to these spin strings the standard approach assumes "fc — A strings" , which should describe the charge 
bound statesd. In the simplest case one assumes two complex momenta k^ are grouped together with a certain spin 
rapidity A from the set of the A 7 : 

sinfc ± = ATij + 0(e-^ L ), (7) 

which one proceeds to insert into the Lieb - Wu equations. This procedure is assumed to hold for k — A strings of 
any length. 

We shall argue in what follows that this approach is flawed (a brief account was given in [ 0) and will introduce 
in its place the BAC approach, which does not rely on a string- hypothesis. We shall argue later that solutions of 
type (0) do not exist for finite L. This is due to the fact that the elementary bound states of (0), to be constructed 
below, cannot be defined with periodic boundary conditions in the presence of unbound electrons. Therefore, such 
states - expected to exist in the infinite volume limit - cannot be obtained from finite volume considerations. To 
obtain a consistent solution at finite volume we need to treat the bound states on an equal footing with the plane 
wave (scattering) states in the Bethe ansatz, rather than trying to recover them as special solutions, "fc — A strings", 
of the Lieb-Wu equations, which were derived from a Bethe ansatz based on plane waves. 

It is amusing to note that already the expression (]?]) hints that it should not be treated, as is conventionally done, 
on par with the spin strings. In contrast to the latter the k — A strings are driven to their asymptotic form in the 
limit L — y oo, i.e. in the infinite volume limit without the need for N to be large as well. This suggests that they are 
not due to a many-body correlation effect as the spin strings (^) and should therefore exist already in the two-particle 
sector of the Hilbert space. 

We proceed now to explain our approach in detail. 

A. Elementary bound states and their boundary conditions 

The Schroedinger equation of the Hubbard model for N particles reads, 

N 

y^Fgx-aN ■■ ■ ) n i ~ 1) •• -, n N) + ■£'ai...(w( Jl l> ■ • • , n, + 1, . . . , U N ) (8) 

i 

N 

+U >J <W.j F ai ... aN {n 1 , . . .,n N ) = EF ai ... aN (n 1 , . . .,n N ). 

i<j 

The solution of (g) in the two particle sector consists of: 

(i) combination of plane waves F aia2 (m, n 2 ) = A aia2 e l ^ kini+k2n2 " > describing unbound particles. These states exist on 
the infinite line as well as on a finite ring of length L, 

(ii) bound state solutions which take the following form on the infinite line: 
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F h (n 1 ,n 2 ) = A s aia2 e^ ni+n ^e-^ ni ~ n ^. (9) 
Here A s aia2 denotes a spin singlet, i.e. A s aia2 — —A s a2ai . The parameters q and £ are related by, 

11 

sinh£(g) = . (10) 

4 cos q 



The energy of the bound state is E(q) = — At cos q cosh £ and £ > 0. The last condition together with ( jl0| ) sets a range 
for g in the interval [— 7r,7r]: 

C/<0: M<f. [ ' 

The signs of the energy depends directly on U: it is negative in the attractive case (U < 0) and positive in the repulsive 
case. In the former case the real momentum q (which is half the quantum mechanical (or crystal) momentum: p = 2q) 
lies near the center of the Brillouin zone and in the latter near the edges. Moreover there is a gap in the spectrum 
E(q) = C/cotanh£ (q): 

U>0: E(q)>\U\ 

U<0: E(q)<-\U\. [LZ) 

The limiting values i\U\ correspond to q — > ±(w/2) from above or below according to the sign of U . In both cases 
£(g) tends to infinity at these points. The state with q = ir/2 (equivalent with q = — tt/2) is strictly local - with a zero 
width wave function both electrons being one on top of the other. The wave function of the local pair is characterized 
by the crystal momentum p = it and energy E = U. This local bound state is the only eigenstate of ([!]), which exists 
also on a finite ring with even number of lattice sites and plays an important role for the SU{2) c h a rge symmetry, see 
below. 

The bound-state solution (||) corresponds, as it should, to a pole of the appropriate S-matrix (or a zero in its 
inverse). We proceed to make it manifest. 

The scattering matrix between two unbound electrons, 

sinfe; — sinfcj + i^Pij 

Sij(ki, kj) = : ; ; ; ■ — , (13) 

smfcj — smkj +i\ 
is used to construct the conventional scattering solution, 

F(m,n 2 ) = Ae^ kini+k2n2 \A aia2 Q(n2 - m) + [S 12 A] aia2 Q(ni - n a )) (14) 

where A denotes the antisymmetrizer and the momenta are real. 
To cast expression @ into this form, rewrite it as, 

A s aia2 e l{q ~ li)ni e l(q+li)n2 0(n 2 - m) + i^/'^^'e^'"^^ - n 2 ), (15) 



which is of the form (14) if we identify: k\ — k~ = q — i£ and k 2 = k + = q + i£, provided [S , i2^4 s ]aia 2 = 0- Using now 
the bound state condition (|lC)|), we find, sin/c 1 * 1 = sing cosh £ ± icosgsinh£ = sing cosh £ ^fij = 4>{q) T ij ■ We have 
used the following notation: 

/ u 2 \ 1/2 

0(g) = K(sinfc ± ) = sing 1 + — 5- • (16) 



16 cos 2 g 



Thus, inserting these values into eq.(|i3|) we have, 



S 12 (k-,k+) = ^(l + P 12 ). (17) 
Observe that S« vanishes indeed for the singlet (P\ 2 = — 1): [*S l i2^4 s ]aia 2 = 

0, while, S 21 = S^ 1 = S 12 {k+,kr), is m 
turn undefinedt^l: the complex momentum satisfying the bound state condition ( |l0| ) is placed on the pole - a standard 
result of scattering theory. 

Hence ( |i4| ) gives correctly the wavefunction (|9|) after antisymmetrization as the amplitude e^ x which diverges in 
interval x > is projected out: 
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F(m,n 2 ) =Ae i ^- i ^ ni e^ q+ ^ n2 A s aia2 Q{n2-ni) = A s aia2 e iq(ni+n ^er^ ni ~ n ^. (18) 

We conclude that the bound state can be written in the usual Bethe ansatz form ( |l4| ) if we allow for singular S-matrices 
like ©. 

If we seek an analogue of this bound state solution in the framework of PBC we find that the momentum is not 
placed at the pole of the S-matrix. Assume two electrons on the ring with L sites. Let us consider a singlet wave 
function in the two-particle sector, N = 2, and 5 = 0, (i.e. M = 1). The equations (0,0) read, 



e 



A — sinfcj + ij 



h J = 1,2 (19) 



2 

A — sin ki — It , , 

1 1 (20) 

A - sinfc,- + if ' 



Now we look for solutions of ( J19| , |20| ) with complex kj, of the form fc 12 = g ± Inserting into the equations we find 
that the spinon parameter A is given by A = 4>(q) and £ satisfies the following equation 



-2££ /cos<7sinh£ 



v cos q sinh £ — | 
which means, 



(21) 



cosq| sinh£ - ^ - e^ L (22) 

for finite L. It follows that the bound state momentum in a finite system with PBC according to the Lieb-Wu 
equations deviates from the pole in the corresponding S-matrix by a term of order e~^ L . This makes S12 regular and 
leads to a F{n\,n2) periodic in each variable separately, because in each sector of configuration space there is also 
an amplitude with the exponentially diverging solution. These amplitudes are forbidden only for L — > 00. The fact, 
that the bound state parameter does not satisfy the pole condition ( |To| ) for finite L renders the bound state unstable 
when scattered on additional particlestfl To avoid this problem, we have to introduce boundary conditions, which 
allow the state (O) to exist in a finite system. This is done in the next section. 



B. The Composite Boundary Conditions (CBC) and the Bethe Ansatz for Composites (BAC) 

We proceed to introduce boundary conditions which allow us to incorporate bound-states consistently, bypassing 
the unsatisfactory "string hypothesis" . We shall show that the Composite Boundary Conditions (i) define a complete 
Hilbert space spanned by 4 L states for finite volume L, (ii) lead in in the infinite volume limit to the conventional 
quantization of the infinite line, (iii) respect the formation of the bound states. We conjecture that they allow, already 
for finite volume, a symmetry expected to hold only in the infinite volume limit. The presence of this symmetry may 
underlie the exact form of our Bethe Ansatz equations. 

To impose a boundary condition on particle j we have to find a path in configuration space leading from region 
Xj < x\ < X2 ■ ■ ■ < xn to region x\ < x-i ... < xn < Xj. If j belongs to a bound state, the corresponding product of 
S-matrices is necessarily singular. We can avoid the singularity by taking the two members j~ and j + simultaneously 
around the ring of circumference L, by imposing the following modified boundary condition for bound states, 

F{x\ . . . , Xj- , . . . , xj+ , . . . xn) — F(x\ . . . , Xj- + L, . . , , Xj+ + L, . . . xn) (23) 

corresponding to periodicity of the center of mass coordinate. For unbound particles we retain conventional periodicity. 
Two new nonsingular S-matrices arise. The S-matrix of a bound pair with an unbound particle, S^._ and the 

S-matrix, between two bound pairs. These are in addition to the usual Sfj 1 given by jl^ ) valid for both 

ki and kj real. 

The factorization property of the S-matrices (|l3|) makes the following construction possible: Because of the validity 
of the Yang Baxter relation it is sufficient to consider only amplitudes in which the two members of the bound pair, 
with momenta k = q ± i£, are neighbors. Obviously the scattering matrix of the pair off an unbound particle with 
momentum k is, 
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^i(23) — $13 (k, k + )Si2(k, k ), (24) 



1(23) 

and we find (see appendix A), 

, sinfc — 6(q) — if 

*»><*■«> - azr^+if - (25) 

S ,ub is a scalar in spin-space, and acts as a pure (momentum dependent) phase shift on the wavefunctions. The bound 
state is stable, its internal wavefunction in is not affected by the scattering with the unbound particle. In other words, 
it couples directly only to the charge degrees of freedom via the phase-shift (p5[). In an analogous way we derive 
S {12) (34) (Qa,Qb) by using, 

^(12)(34) = ^1(34) 

(K,q b )S^ 34) (k+,q b ) (26) 

with the result 

S mm (q a ,<l b ) - ^ qa) _^ qb)+q - (27) 

We proceed now to derive the Bethe-Ansatz equations. As S ub and S bb commute with each other and all the S uu, s 
they do not enter the self consistency BA equations (^). This fact is responsible for the (partial) decoupling of the 
bound states from the free particles which are correlated via the spinon parameters A~. Assume N = N u + 2N b 
electrons, where 2N b particles are in bound states characterized by momenta qi, I — 1 . . . N b and N u particles are 
in plane wave states given by momenta fcj, j = 1 . . . N u . The total spin S is given by S = h(N u — 2M), where M 
denotes the number of spin-lowering operators in the algebraic Bethe ansatz. 

One proceeds to construct the Z- matrix which takes a particle or a bound pair around the ring L. For a particle j 
in an unbound state it takes the form, 

ry QUU QUU QU b QU b /no\ 

* 3 — D lj ■ ■ ■ o N ujOij . . . J N bj (40) 

with eigenvalue Zj = e lkjL . The last N b S-matrices are phases. Diagonalizing Zj is a standard procedure and we find 
the BA equations, 



yI Xs — sin kj — if t~t 4>{qi) — sin k 



tt ^ aillf "J ~ 4 TT VWJ Z 1 ~ 4 (29) 

s=i Xs ~ sin fc j + * f r=i tin ) - sin k J + i i 



-pj- A 7 - Xs - 1 5 = -pr A 7 -sinfc J -t I 

lA Xi _x 5 +q 11 A 7 - sinfc, + if [ > 

Similarly the Z-matrix for a bound pair determines its momentum, and is given by, 

y Qub Qub obb Qbb /oi \ 

- b m) ■ ■ ■ b N^ij)t> {r - <r + ){ij) ■ ■ ■ *(r- b ,r+ b m ^ 

with eigenvalue Zuj) = e 2lqiL . The resulting BA equations read, 



n X sinfcj- -4>{qi) -if -pr <j>(g n ) - (j>{gi) - if 



=i sin k 3 - (j>(q{) + if 11 4>{q n ) - (f>{qi) + if 



The equations ( |29| , |30| , |32| ) replace the set (|3|,|4|) in the presence of two-particle bound states. These bound states, 
however, are not the only multi-particle complexes which are allowed by the hamiltonian. We can infer the existence 
of higher composites by looking for zeros and poles of the S-matrix for two two-particle bound states (E7|) . The zero 
of the bound - state S-matrix, (p7j), is at, <f>(q a ) — <fi{q b ) = if, corresponding to complex g's. Choosing <p\ 2 = 0o ^f , 
we find the four momenta of this "double" bound state, the so-called quartetB, 



G 



sinkf = 


<t>l 


.u 
+ 1- 

4 


sinfc| = 




. u 

— i — 
4 


sinfcg = 




u 

+ *4 


sin k\ = c 


t>i- 


u 



(33) 



In other words, 



k{ = 7r — arcsin(0 o + i — ) 

k\ = arcsin^ 2 ,) 
fc 3 = 7r — arcsin^ 2 ,) 

k\ = 7r — arcsin(0Q — i—). (34) 

Again we can derive the S-matrix of this state with the unbound and the simple bound state. The latter follows from 
the former upon application of a relation similar to ([26]) . We find 



w _ sinfc-<ftg-if 



2 



In general there are bound complexes of 2m electrons, (m-complexes) , corresponding to a pole in the S-matrix of a 
(to — l)-complex with a simple bound pair, (a 1-complex). They are parameterized by to complex numbers {</>™, 
7 = 1... to} in the complex plane: 



corresponding to 2m complex momenta, 



+ (m + l-2j)i- (36) 



u 

k™ = 7r - arcsin(</>™ + im-) 

11 

k? = arcsin(^ + i(m - 2)-) 

11 

fc™ = 7r - arcsin(0™ + «'( m _ 2 ) ^) 



(37) 



v 2m-l 



7r — arcsin(0Q — i(m — 2)—) 



fc 2m = 7r - arcsin^ 2 , - i(m)-) 
In fact, these correspond to the bound complexes conjectured by Takahashii. 

We proceed along the same lines as before: The S-matrix of an m-complex with an unbound particle is, 



ju(m) _ sin fc - <ffi - imf 
^ sin k - 0™ + im f 



Q u(m) _ -""4 ,„„n 

^fc^ - 77T7 ISTmr 



allowing us to derive the S-matrix of an m-complex with an n-complex: 

g (m)(n) _ <K~ M - I" - m\q ^ - -jn + m)q min( fi n) ~* / flT-jg- (jn - m| + 2Qg \ 2 

^Oi^o* A™ _ At? -I- In — m\i— rh™ — rh 7 } 4- (<n -i- m\i — H \ MP> - rh 7 } 4- (\<n - m I 4- MUM I ' ^ ' 



+ | n _ m |i» 0™ - 0" + (n + m)i| ^ _ + (| n _ m | + 2 /)i» 

If we impose the composite boundary condition on the TO-complexes, 

F(x 1 ...,xt\...,x i Z ) ,--^N) = F(x 1 ...,x ( r ) +L,...,x^+L,...x N ) (40) 
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where the {x^} are the coordinates of the members of the m-complex, we are led immediately to the corresponding 
eigenvalue of the m-complex transfer matrix, (recall, all bound states are spin singlets): 

exp(iL^ k^) = exp(iL[-29?arcsm(>o + rra~) + (m + 1)tt]) = e iqm{ML . (41) 
i 

In this way we get the full BAC equations, a generalization of ( ]29| , |30| , |32| ) : 

ik } L _ TT A 5 -sinfc 3 - i| „ ( )u 

TT ^7 ~ - if _ T-r A 7 - Smfcj - I j 

N" 

k i<P 



e ^i^ = H S^I^US^ (44) 

(0,71)7^(6,771) j 

where the index (a, n) runs over the set of all n-complexes present. These equations were also derived by Takahashi 
in [ ||, within the framework the k — A string hypothesis. His procedure, however, involves discarding finite volume 
corrections in equations written for finite volume. Our derivation makes no use of this hypothesis as all charge bound 
states are consistently incorporated ab initio. 

To discuss the nature of the cigcnstates and count them it is convenient to consider (0) in a logarithmic form, 



b { b ' m) )L 



£ e nm (^- m) -4 a - n) )-^J b m (45) 



where, 



e nm ( x ) = e( r ^-)+26(- — x -— +...20 (- — x -— )+0 (-£-) „ ,n 

\n — m\J \\n — m\ + 2 J \ \n — m\ — I ) \n + m 



*(f)+2*(f) + ...2* 



X 



e(— 

2n-2) ' \2n 



with 6(x) = -2tan" 1 (x^). 

Each .allowed choice of quantum numbers {J^} uniquely labels the eigenstate, and the allowed ranges can be 
deduced^ from eqn(^5|), (a derivation for the case m — 1 is given in the next section), 

| J* | < hL - N u - ]T tnmMn)- (46) 

Here M n denote the number of n-complexes, t nm is defined as t nm = 2min(n, m) — 5 nm . These equations are the 
starting point for counting the number of states of the model. 

We end this subsection by showing that the dimensions of the Hilbert spaces for CBC and PBC are the same. The 
dimension of the Hilbert space with CBC is larger or equal to the dimension of the space with PBC: 

dim TLcbc > dim TLpbc (47) 

as each vector in TLpbc lies also in TLcbc- On the other hand, there is a one-to-one relation between states in TLcbc-, 
not satisfying PBC and states in TLpbc- If one writes a vector in TLcbc as function of the center-of-mass and relative 
coordinates 

F(xx, ■ x N ) = F(X c , m ,,Xi ~ Xj) 

one has the injective mapping &(F) onto TLpbc'- 

$(F)=F(X e . m .,[x i -x j ] modi) (48) 

because the value of the wavefunction F(xi, . . . , xn) outside the interval \xi — Xj \ < L/2 is completely determined by 
its value inside. It follows 

dim TLcbc — dim TLpbc (49) 
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C. Global Symmetries and the Bethe-Ansatz nature of all Eigenstates 



We begin by considering the global symmetries of the model whose eigenvalues (partially) label the eigenstates. 
We shall show, among other things, that all the eigenstates of the model are of the Bethe Ansatz form. This is in 
contrast to claims by Efiler, Korepin and Schoutensa that only the highest weight states have the Bethe Ansatz form 
while the rest cannot be represented this way. 

First, we have the spin S'f7(2)-symmetry, with 

i i i 

The spin operators commute with the particle number operator, N = + n i,ii ana - ^he application of the 

spin-lowering operator S~ does not lead to a change in particle number. 

The spin highest weight state \S = S z = (1/2) (JV - 2M)) is defined by a solution {A 7 , 7 = 1...M}, all A 7 
finite. To obtain another member of the multiplet consider the solution with M + 1 spin rapidities {Ai . . . \m, Ao} 
and Ao = 00. It formally satisfies ( |4^^ ), and corresponds to the state obtained from {Ai . . . Aa/} by an application 
of S~ , since -B(A) — ► S~ /A, when A — > 00, in the framework of the algebraic Bethe ansatz. Repeated application of 
S~ will generate the whole SU(2) spin multiplet from {A 7 }, which consists of 25 + 1 states, i.e. we will find a zero 
wavefunction for more than 2S spin rapidities equal to infinity. A consistent way of defining this process is provided 
by considering an anisotropic version of the spin sector. New complex roots, denote them as A - - roots, are generated 
with 7r/7 as their imaginary part, where 7 is the anisotropy parameter. When 7^0 isotropy is regained and t h^A-^ 



An even simple^ argument holds in the case of the charge SU (2) symmetry. It is defined on a ring with L sites (L 
even), as follows! 



roots are sent to infinity in a controllable way, generating the Bethe Ansatz states that complete the multiplet 

ti, = \{L-N), r,+ =^(-1)%^, ^=E(- 1 )^iT' ( 5 °) 

3 j 

The algebra is consistent with the CBC since the 77 ± operators create (destroy) local pairs - see below. The 
symmetry is manifest when one adds a chemical potential term with fi = — U/2 to the hamiltonian: 

H h . f .= -^L+iVw + h.c.) + Un hl n Lt - —N (51) 

1— — OO 

This choice of \x corresponds to a half-filled system in the grand canonical formalism. We have 

[H h . f .,tf k ]=0 [JV',fr t ]=±2» 7 ± (52) 

Clearly, the symmetry generators rf^ 1 mix sectors with different particle number. 

In terms of the BAC equations the ^-symmetry has a simple explanation: Consider a given eigenstate of H , say 
hA) = IV'({A 7 }, {kj}, {(#})) with N particles characterized by unbound momenta {kj} and bound momenta {qi}, as 
well as spin content given by {A 7 }. Acting with the operator r] + adds to it a bound pair with q = n/2 and crystal 
momentum p = tt: 

r,+ m{X 7 },{k,}AQi})) = W({\yh{kj},{m,n/2})}. (53) 

(We have assumed qi ^ 7r/2, VZ). Note that for a finite system, namely with CBC for L < 00, the state with q — tt/2 
exists for even L. The state |^') is again an eigenstate of H, because the S-matrix between the bound pair with 
q = n/2 and all other complexes is the identity (total transmission), - see (|38|), ( ]39| ) and recall: <P(tt/2) = 00. The 
state \ip') has then N' = N + 2 particles and its energy is 

H\iP') = [E{ip) + U]\i)') (54) 

since adding a bound pair at the edge, q = tt/2, corresponds to AE = U, see eqn ([l2|) and subsequent discussion. 
Thus, 

H h4 .W) = E h .f.{^)W)- (55) 
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is degenerate with \ip) and the symmetry is manifest. We now proceed to add several pairs with q = ir/2 to the 
state \ip). These states have zero width - £(±7r/2) = oo - and they behave as hard-core bosons: The S-matrix among 
themselves is S = — 1, corresponding to total reflection. The maximal number of applications of rj + to \ifi) is restricted 
to L — N, the number of available lattice sites. Here N includes bound as well as unbound states - the former with 
q ±7r/2 and thus £ oo, giving them a finite spread. For more than L — N applications of 77 + the state is annihilated 
due to the Pauli principle. The total number of states degenerate under the 77-symmetry is therefore L — N + 1. This 
degeneracy coincides with the dimension of the SU(2) c harge multiplet for S c harge — S z<c harge = \{L — N), which is 
the eigenvalue of r\ z applied to \if>). It follows that the Bethe state |?/;Vis a lowest weight state of this symmetry, and 
all members of the multiplet have the appropriate Bethe Ansatz formEil. 

Actually, our construction goes further. We see from (|54|), that the 77-symmetry can be used to group the eigenstates 
of H into multiplcts even away from half-filling when the SU(2) c ] large symmetry is explicitly broken: The energies of 
the L — N + 1 states in the multiplet are equally spaced with E^+i — Ei = U — 2/i. 



D. Periodic Boundary Conditions (PBC), and Takahashi's conjecture 

In this section we shall discuss the approach to the bound states within the usual scheme - imposing periodic 
boundary conditions. 

Defining the system on a finite ring L and imposing PBC (a multi-torus geometry) we have seen above that the 
parameters q and £ of a normalizable eigenstate of H in the two-particle sector deviate from the pole (given by eq. 
(|l0|)) in the scattering matrix by a term of order e~^ L : 

I cosq| sinh£ - ^ - e~ ?L - e~^ L for \q\ resp. \q - ir\ < n/2. (56) 

Tlus deviation is responsible for the instability of this state when a third particle is added with real momentum 
kM. The S-matrix S ub is no longer a pure phase, and alters the spin structure of the bound state. There is a finite 
(although exponentially small) probability for the bound state to switch from the spin singlet to the spin triplet state 
upon passing through the third particle. As the wavefunction has to be antisymmetric in configuration space, we 
have an anti-bound state which cannot satisfy periodic boundary conditions and is clearly forbidden in the infinite 
volume limit. That means that PBC contradict the local properties of the interaction encoded in the pole structure 
of the S-matrix, while the CBC are consistent with it. The coordinate space of the system with the CBC is not a 
multi-torus, but has a more complicated topological structure which, however, turns into the infinite line for L — > 00 
as does the geometry of PBC. 

One might still look for the analogues of the normalizable two-particle state (and higher composites) on a finite 
ring with PBC - solutions to the Lieb-Wu equations (||,||) with complex momenta. It is clear that some of the states 
must contain complex rnamenta because otherwise the subspace containing doubly occupied sites would be projected 
out in the U — > 00 limitEJ. n 

Takahashi's k — A string hypothesise! assumes that the states containing complex momenta are of the special form 
for large L, i.e. they become elementary bound states for infinite volume. Therefore, using the string hypothesis 
one obtains BA equations similar to those for CBC (^2|,fl3|,fl4]) , but containing correction terms. In fact, if one drops 

the corrections ~ e~~ L then Takahashi's equations do not describe PBC but instead CBC for finite L. The proof of 
completeness of the Bethe ansatz solutions given in [ uses Takahashi's equations for finite L and counts not the 
number of states on a finite ring for which they are approximate but on the CBC geometry, where these equations 
become exact. But is Takahashi's hypothesis correct for PBC and finite L when the terms proportional to e~~ L do 
not vanish? 

We argue now, by examining the consequences of this assumption, that this is not the case and the complete 
spectrum for PBC is (in general) not given by k — A string solutions to (|^,^) . 

We begin by reviewing Takahashi's approach in a simple case with no spin-strings or higher composites present. We 
assume N e electrons 2N h of them carrying complex momenta = qi ± . In addition we assume M — N b + Al u 
real spin rapidities {A;, A 7 ; I — 1 . . . N b ; 7 = 1 . . . M u }, N b of them, {A;}, associated with the complex momenta 
as follows, 

sinfc± ee 4>{ qi ) =f 1x1 = h T i\ + 0{e-^ L ). (57) 
Plugging this form into (|],[|) one finds, 
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11 X _ sin £• . 4-i™ 11 A , — sin 4- i2i { -°° ) 



sin /c., + «t A; — sin fco 

J 4 i =1 1 J 



A7 Aa 2 7^ -| — r A7 sin kj i^j 



tt- A; — Aa — ^-r A/ — sin A;, — i% . , 

nfeit-n ,,-^^ ^) <»> 



A 7 - A; + if 4- A n - A; + i| 

7=1 ' z ri^Z z 

Here {£} denote terms of order e _ftL , with n > |||. The phase e"^ ! is defined as 

A; — smkt — i¥ 



A/ — sin fc ; 



(62) 



Now the set (||8| - ^l|) consists of N u + M u + 2N b coupled algebraic equations for the variables {kj, A 7 , qi, A;}, i.e. 
the number of variables coincides with the number of equations^. But in the L — > oo limit, A; is not independent 
from qi, because of (f>7\): The A; should be eliminated from the set together with the N b parameters ipi, which do 
not describe physical properties of the state. Takahashi did this by substituting eq. ( |60| ) into ( |6l| ) after dropping the 
£-terms. The resulting equations are the set ( |29|j30| , p2] ) . It is clearly consistent, as it can be derived using BAC. The 
question arises whether the finite size correction terms in (|5^ - |6~i| ) can spoil the consistency for finite L. We proceed 
now to show that this is the case in general. 

Let us assume there is a consistent solution {kj(L), A 7 (L), qi(L), A/(L), ifji(L)} of (^8|- ^l|) for arbitrary (large) L. 
We define the set {fc°(L), q?(L), A° (L), A?(L) = <j>(q?(L)),tpf(L)} as the solution of the zcroth order (in £) terms of 

(§-13): 



^ A°-sinfc°-zf iL 0( g °)- S infc°-zf 
11 A o - sin fc° + q 11 0(g?) - sin fco + z| 



TT A 7 ~ A a ~ 1 2 _ TT A 7 Sln S" Z 4 

Al AO -Ag-h.^ HAO-sinfeO + il l ^ 

11 0(g o) - A° + if 1 = 1 - sin *° + if 6 l ° 0j 

e ^o L = fr A 7 -^°)-^ ff #9°)-^)-^ , 

11 A, - *(g?) + if 11 «%°) - Hq 0) + q ■ W 

The set ( |63| - |66] ) contains again the same number of unknowns as equations and should give the same solutions 
as ( p9| , |30| , [32| ) for the parameters {A°}, {q^}- In addition it determines the phases e 1 ^ 1 . Because we dropped 

only the exponentially small correction terms £, the set {kj, A°, qf, A = 4>(q^)} deviates from {kj, A 7 , qi, A;} only in 
quantities of order £ . Especially, if we choose an arbitrary I and define e; = e~^ lL , we can write 

k j {L) = k°j{L)+kf\L)e l 
\ g (L)=\° g (L)+\W(L)e l 

q l (L)=qf(L) + qi 1) (L)e l (67) 
A l (L)=<f>(q°(L))+AW(L)e l 

Xi (L) = ^ + x ; (1) (^. 

We will now show that the coefficient Aj in (^) is determined through two independent equations, leading to an 
over-constraint. 
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One set of equations to determine is obtained as follows: dividing (||) for with (|^) for k x we get the following 
equation: 



-2(iL _ 



(A, - 0(gQ) 2 + (6 - | ) 2 (A 7 - 0(g,)) 2 + (ft - \f 1f (A» - 0( g; )) 2 + (fr - |) 2 

(A, - 0(<zO) 2 + (6 + f ) 2 (A 7 - <t>( qi )Y + + f ) 2 W (A„ - 0(*)) 2 + (6 + f ) 2 ' 1 J 



which leads to, 




s #i _ i 
2 iV>? + 1 



(x, (1) ) 2 /({A 7 ,g°}) 



(69) 



with / some function of {A 7 , q^}- This equation determines x ; as function of the zeroth order variables {A 7 , q^, itf}, 

therefore also the coefficient Aj in (|67j), (using (p2|)), up to terms which are by themselves exponentially small. On 

the other hand, we obtain A^ directly from eq. (|60|), by looking at the ej-correction. These two determinations of 

the coefficient A^ are independent, as the former is derived from (0) and the latter from (^). This indicates that the 
set (p8| - |6l| ) is over-constrained in the first order finite volume correction to the thermodynamic limit. This effect can 
be studied in the three particle case by an explicit construction of the BA wave function assuming the k — A string 
hypothesis. An over-determination of the parameters in expressions of order £ is foundta. We conclude that the fc — A 
string hypothesis does not correspond to an actual solution of the BA equations for sufficiently large but finite L. The 
spectrum of the model on a finite ring is not in analytic one-to-one correspondence with the spectrum on the infinite 
line. On the other hand, the CBC - spectrum develops smoothly into the infinite volume limit, probably because it 
already contains an additional symmetry. This symmetry is presumably destroyed by the PBC and appears in this 
case only in the infinite volume limit. 

Another way to observe the over-determination of (|],f|) for finite volume is based on the algebraic Bethe ansatz. 

We show that equations (|60|), which do not contain L explicitly, become redundant for L — > oo, leading to superfluous 
constraints on the parameters: assume a pair of complex conjugated momenta k + , k~ related to the spin momentum A 
by (0). An eigenstate of the corresponding inhomogeneous transfer matrix Z(fi) with (arbitrary) spectral parameter 

ix then has the form (see e.g. [ |): = FJ^i B(X 7 )B(A)\ui) for total spin S = ±{N U - 2M U ), where the M u 
creation operators 5(A 7 ), acting on the ferromagnetic vacuum create M u |-spins, and B(A) creates the J,-spin 
of the bound state. By explicit calculation we find that B(A) acting on \oj) diverges as e~^ L for L — > oo, whereas 
the .B's not associated with the complex pair do not diverge. Normalizing \ip) by multiplication with e~^~ L yields 
accordingly an exponential suppression of vectors of the form Y[s B(\s)\uj) which do not contain B(A) among the 
B(Xg). Now the equation ([30]) for A is necessary to cancel an "unwanted term" in the eigenvalue equation for 



(Z(li) - E(p)M) 



7 



M" 



M" 



JJ B(X s )B(A)B( f i)\Lu) + a Q JJ B(Xy)B(fx)\u;}. 



(70) 



5^7 



It ensures, in particular, that ao = 0. But, as argued previously, the vector FJ B(X 7 )B(fj,)\uj} is projected out in 
the infinite volume limit (exponentially suppressed for finite L), and therefore ( J60| ) is not necessary for to be an 
eigenvector of Z(ll) for L — > oo. That means that only a subset of all states allowed in the infinite volume can be 
generated by the L — > oo limit of string solutions to (|3|,|]) for finite L. 

In conclusion we find that the k — A string hypothesis represents an over-constrained ansatz for solutions of BA 
equations for periodic boundary conditions, having in general no solutions for large but finite volume L. 



III. GAPPED EXCITATIONS 



In this section we will compute the simplest gapped excitations above the ground state for the repulsive and the 
attractive case. In the former case this amounts to the formation of a bound pair above the sea of unbound particles, 
in the latter case one of the bound pairs forming the ground state is broken and the resulting excitation has nonzero 
spin. Both excitations are characterized by a gap of order U and are accompanied, when the number of electrons is 
held fixed, by two gapless excitations: two holons in repulsive or two dressed electrons in the attractive case. 

The range of the gapped excitation momentum, p hs , depends on the filling n, — 7r(l — n) < p bs < ir(l — n). At half 
filling it becomes therefore a non dynamic mode. We present the dispersion of the gapped mode for various fillings 
and interaction strengths. 
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A. The repulsive case 



A single bound pair above the ground state of the repulsive Hubbard model can be created without changing the 
particle number N by placing two holes in the sea of the charge quantum numbers. N u = N — 2 and M = y — 1 in 
the notation given above. Equations (^,[^,[32]) then reacO, 

M 

Lk-j = 2nrij + 9i(sin kj - Xs) + 9i(sin kj - <j>(q)) (71) 

6=1 

N-2 M 

^2 Si(sinfcj -Xy)=^2 Q 2 (Xs - A 7 ) + 2ttI 7 (72) 
3=1 s=i 

N-2 

2qL=^2e 1 ((f>(q)-sink j ) + 2TrJ (73) 

3=1 

where Q n (x) = 9(x/n). The range for the quantum numbers {rij} is: —N/2 < rij < N/2 — 1 for M = N/2 — 1 even and 
-(N-l)/2 < rij < (JV-l)/2 if A/ is odd. The {/ 7 } range between —{{N — 2)—M — l)/2 and +((iV — 2) — M — 1)/2. 
The {rij} sequence contains two holes as the actual number of free fc-momenta is N — 2. The {/ 7 } sequence does 
not contain holes in the absence of spin excitations. J, the quantum number associated with the bound state, is an 
integer if N is even and a half-odd integer if N is odd. We assume N even in the following. To find the limiting 
values for J, we consider the boundaries of the allowed range for q: ir/2 < q < it, (resp. — it < q < —tt/2). We may 
treat the range for q as connected by shifting — tt/2 to 37r/2. Setting q — tt/2: 

ttL = 27T,r - (N - 2)tt (74) 

as 4>(tt/2) = oo. For q — 3n/2 we have <j>(q) = —oo. Thus, 

3ttL = 2nJ + + (N - 2)tt. (75) 

Hence, 

N N 

L- Y >J> Y (76) 

At half- filling: (N/2) < J < (N/2), which leads to maximal restriction of the phase space for the bound state. In 
the infinite volume limit which is treated here, the phase space for J vanishes at half-filling. RecalLthat this is the 
case for spin singlet excitations where the string parameter is given in terms of the hole parametersO. The Lieb-Wu 
equations studied in [[ll]] lead to a similar rigid relation between q, k\ and ki, namely <fi(q) = ^(sinfci +sinfc2). While 
this rigidity holds at half-filling it is physically meaningless, as the physical momentum p bs becomes independent of q 
at this point (see below). _ . 

We proceed now in the usual manneiO, using the notation of [|||. One introduces the functions {p(k),a(X}} to 
describe the densities of the {k} and {A} solutions respectively. The ground state densities {po(fc), co(A)} are solutions 
of the appropriate integral equations with no holes in the distribution of the kj , while the state under consideration 
now is given in terms of{pb(k) 1 ab(X)}. The integral equations for these densities read, 

Pb (k) + h(k - fcO) + - fc 2 ) = 



- — h cos k 
2ir 



d\a b (\)Ki(smk - A) + j- cos kKi (sin k - <f>(q)) (77) 



cr fc (A) = J Q 6kp b (k)K x (sink - A) - J <LV 'a b (X')K 2 (X - A') 



(78) 



where fci, fea denote the positions of the holes, and K\, K2 as well as R, KP (see below) arc kernels of integral operators 
defined as in [ H . The parameter Q defines the range of k and is given by, 
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1% ^ = f 
I Q Q dkpQ =v 

it depends therefore implicitly on k%, ki and g. The equation for q reads, 

2tt r Q 

2q=—J + dkp(k)Qi(4>(q) -sinfc). 
It is convenient to introduce p'i(k), ui(A), 



p(fc) - plftfc) + ipi(fc) - S(k - hi) - S(k - k 2 ) 
a(X) =a (A) + |(T 1 (A) 

these correspond to the density changes induced by the excitations with respect to the density 
transform of the spin density ci(A) reads in terms of p'i(fc), 

1 rQ 11 i ii 



*i(p) = o Pi(fc)sech(-p)e-^ sin ' £ - -sech(-p)[e 



2J Q ^ ' 2 M J 

It is possible to split p[ (fc) into three parts 

p[(k) = pUkM) + P§(fc,*a) + P 6s (fc,g) 

such that 



4/4 

/C Q [p!j](fc) = — cos k—R I — (sinfc — sinfc,- 
J u \u 



for j = 1,2 and 

/C°[p 6s ](fc) = cos fci^ (sinfc - 0(g)) 
The excitation energy Ai? = £(Q) — Eq(Qo) is defined via 

£ (Q) = -2tL/ Q Q dfcp^(fc)cosfc 
£J(Q) = -2i£ J Q Q dfcp(fc) cos fc 

With the usual definition of the chemical potential /i 

(n\ ( dE A f dE o(Q) \ f 9N (Q )Y 1 



and 



It follows 



rQ 

N (Q) = L / dkp$ 
J-Q 



/Q rW 
dkp[ cos fc + 2i[cos ki + cos fo] — p(Qo) / dfcpj 
-Q ~ ' * ./ Q 

/i is independent of k\, fc2 and g and given by Coil's formulaB 

KQo) _ cos (Qo) - J Q q dfccosfcp c (fc,Q ) 

The excitation energy is a sum of three parts 
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AE = + E\ + E bs (91) 
Each of them consists of a direct, a backfiow and a ground state contribution: 

rQ rQ 
Ej(kj) = 2tcosfcj — 2t / dkp c (k,kj) cos k - fi / dkp c (k,k 3 ) (92) 



The hole-energy goes to —p for kj — > Q. The bound state energy reads 

/■Q /-0 

E bs (q) = -At cos qcosh(£(q)) -2t dkp bs cos k - p / dfcp bs . (93) 

J-Q J-Q 

We proceed to discuss the energy dispersion. We shall do it for any filling. Consider the half filled case: using the 
identity — At cosgcosh(£(q)) = U + 2t J dfccos 2 kKi(sm k — (j>(q)) we see E bs (q) — U vanishes identically at half-filling 
and the dispersion of the excitation depends only on the hole part apart from the gap U. 

Away from half-filling the bound state contribution becomes an independent excitation with its own dispersion. 
The momentum has contributions from the holes and from the bound state, given by, 

AP = P-P = -pf -p%+ p bs . (94) 



The hole- momenta =r-rij read 



V ) = [ 3 dk Po (k) (95) 



(which leads to the identification of the point kj = ±Q with the charge Fermi-momentum ±kp = ±tt(N/L)), while 
for p bs we find 

rQ 

p bs {q) = 2q- dfcp (A)ei((/)(g)-smfc). (96) 
J-Q 

As q runs over the allowed range ir/2 < q < it (resp. — tt < q < —ir/2) , p bs varies between ir + kp and 2n (resp. 
— 2tt < p bs < — 7r — kp). This corresponds to a symmetric band around zero between —\ir — kp\ and tt — kp. It follows 
that at half-filling p bs = 0, i.e. the phase volume vanishes. The physical meaning of the shrinking to zero of parameter 
range at half filling has a simple interpretation: the bound state has no room to|-oropagate. Although the bound 
state parameter q is fixed in this case, in terms of the hole parameters ki and fe'^J, this has no physical meaning: 
Energy and momentum, the physical parameters of the excitation, are independent of the unphysical parameter q, 
namely E bs =_Z7 and p bs = 0. This redundancy of the parameter q at half- filling has led to some misconceptions in 
the literatureEi 

Away from half-filling, the dispersion of the (always gapped) excitation is given parametrically by ( p3| ) and (Q) for 
fixed hole-momenta k\ and k 2 . We solved (|^,|5|) numerically for t = 1 and two values of the interaction strength, 
U = 0.5 and U = 2. Figures 1 and 2 show the dispersion of the bound state energy E bs (p bs ) for U = 2 and U = 0.5 
respectively and for different values of the density n. Figures 3 and 4 show the dispersion curves for the same values 
of U, but for densities very close to half filling. 

We notice that beyond a critical filling n c (U) the bound state energy dips below U. The critical value depends on 
U and decreases with increasing U. We are currently exploring whether the gap between the top of the holon and 
spinon bands and the bottom of the upper Hubbard band can actually close for some value of U. 
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FIG. 1. The band of the bound state excitation for U = 2, and densities n — 0.44, 0.62, 0.78 and 0.85, respectively. 
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FIG. 2. The band for U = 0.5 and three values for the density, n — 0.53 (long-dashed), n = 0.75 (dashed), and n = 0.96 
(solid line). 



1G 



2.20 



o 2.10 




0.0 

momentum 



0.5 



FIG. 3. The dispersion for U = 2 and n = 0.84 (long-dashed), n = 0.87 (dot-dashed) and n = 0.89 (solid line). The critical 
density n c ~ 0.86. 
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FIG. 4. The dispersion for U = 0.5 and three densities close to half filling: n = 0.97, 0.98 and 0.99. The critical density 
n c = 0.98. 
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B. The attractive case 



In the attractive case U < and we expect the ground state to consist solely of bound pairs (for even N), being a 
total spin singlet with N b = N/2 and N u = M u = 0. We find accordingly, 



JV" 



2 lgi _ TT v\Qn) - 4>(qi) - ij 

a reduced form of (|32|). We proceed in analogy with the treatment in section III A. The set of the qi lies in the interval 
-f < qi < f . Taking the logarithm of @ 

N b 

2q l L = 2-Kj l -Y,®2{<t>{<li)-<t>(<ln)) (98) 

The range for the quantum number J; , corresponding to the range for qi above, is 

- (L - N b - l)/2 < Ji < (L - N b - l)/2 (99) 

i.e. Ji is integer (half-odd-integer) if L — N b is odd (even). At half-filling N = L, the J/'s are filling all the slots 
allowed by fl99|). It is convenient to change variables from qi to 4>(qi), in defining the density, La(4>) — dJ/d<fi(q) which 
leads to the ground state integral equation, 

C B " [a] (</>) = -F{4>) (100) 

with L B °[a]{4>) = a{4>) + / S ^ n cl^^ (0 - 0') CT (</>') and tlic inhomogeneous term, F(4>) = dq(0)/d0 = 
- {<f> + ij) 2 }- As observed in [ |, operator C B plays the same role as KP in the repulsive case. The in- 
tegration limit Bq is determined by J"_J cr('/ ) ) = N b /L. At half-filling the r.h.s. is 1/2, allowing us to deduce from 
that Bq — oo, and the equation can be solved via Fourier transformation!. 

We now consider the simplest gapped excitations above this ground state. It involves pair-breaking, and is therefore 
a spin excitation. 

We consider first the triplet: it is created by removing one bound pair and adding two free particles in a triplet 
spin state. We have, 

N b -1 



2q t L = 2nJi - Q 2(^ ~ $ «) - ©i(0J - sinfci) - 0i (0j - sinfc 2 ) (101) 

N b -1 

k 3 L = - ^2 ©i(sinfcj - 4>i) (102) 



for j = 1, 2. The total spin is S = 1, N b -> iV 6 - 1, A^" = 2, M u = 0. 

At half- filling the number of available slots is reduced by one (see (|9^)), therefore no hole opens in the J- sequence. 
The parameter of the excitation are just the two momenta of the free particles k\ and hi. Away from half- filling 
N = 2N b < L and not all allowed slots are occupied in the ground state. Now, Bq < oo, the Jz's are distributed 
symmetrically around zero and |q mQa; | < n/2. That means, we can create a hole in the J- sequence at position Jh, 
which corresponds to a \4>h\ < Bq although the number of allowed slots decreases by one. This Jh is the third parameter 
of the spin-triplet excitation and corresponds to the bound state parameter in the repulsive case. We proceed by 
introducing the densities, a((f>) = a B + j^a\{(j>) and ai((f>) = p'\{4>) —&{4> — 4>h) where a B {(j)), being the ground state 

distribution with Fermi level B rather than £?°,is determined by the normalization condition f_„ cr(4>) — (N b — 1)/L. 
Expressing the smooth density eri(</>) as a sum of three terms, a[(4>) = cr^(0) + cr?(</>) + <Jh{<t>) we find, 

C B [(J 3 C }(4>) = -K x {<t> - sin*,-) (103) 

for j = 1,2 and 

C B [a h ((f>, <}> h )]{dp) = K 2 (4> - <M. (104) 
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Having solved ( 103| ,104) for the densities cr|, er/j, we compute the excitation energy. The total excitation energy consists 
of three terms, 

2 B B 

AE = + E 2 +E h = -2tJ2 cos k J -£(cj) h )+ f ° £((/>) a^) - \B Q ) f ° a'M, (105) 

j—l J — Bo J Be, 

associated with the two unbound electrons and the independent hole respectively, 

Ej = -2tcoskj + [ <r 3 c ((j))[£(<f>) - n a ] (106) 

J-B 

E h = -£{^ h ) + f ° a h {& 4>h)[£{4>) - H°]. (107) 

J -B 

fi a is the chemical potential and accounts for the shift of the Fermi momentum q(B) through the excitation. It is 
defined with respect to the number of bound pairs, [f = dE /dN b , and given by an analog to (|9C|): 



_ g(g o )-/_V(0K(6Bb) 

A 1 (Bo) " jWo — — • (108) 

1 ~ J-B o a h(<P,B ) 



£(4>) denotes the bound state energy function, 



£(</>) = -4^1 -(0 + z J) 2 <-\U\. (109) 
At half-filling <Th(4>, oo) = and /i CT (oo) = — \U\, as expectedi. 

The momenta pf, ph, associated with the dressed electrons j = 1,2 and the hole respectively, are given by, 

p h = f h ^°{4>) (no) 

Jo 

/Bo 
d^6i(sin% - 4>). (Ill) 
-Bo 

Removing the bound state at the Fermi level B we have, 

p h (B)=Tr^ = ~n (112) 

which identifies the Fermi momentum in the attractive Note also that the dressed momenta pj of 

the two unbound elec trons devia te from their free values kj . The hole contribution to energy and momentum, given 
by Eh and ph in (107) and ( |11C| ), play the same role as E bs and p bs in section III A. 

Now consider the singlet excitation. We break a pair and put the two electrons into a spin singlet state. Again N b 
is reduced by one, leading to an additional degree of freedom away from half-filling. N u = 2, but M u — 1 and S = 0, 



the total spin is not changed. We have one A parameter for the unbound electrons. Equation (|30|) reads then, 

^ _ A - sin fci + ijf^ A - sin fc 2 + 
A — sin fci — i^f A — sin k 2 — i^-j- 

which leads to the familiar form 

A = -(sin fci + sinfc 2 ). (114) 
Eq. ( |101| ) remains valid in the singlet case while eq. ( |102[ ) becomes 

Af b -1 

kjL = 2-KUj — 8i(sinfcj — <fii) — 0i(sinfcj — A) (115) 
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The effect of this modification is to change the relative phase shift of the particles, 

gsinglet = gtriplet _ Q^fa ^ „ sin fc 2 ) ( 116 ) 

as expected on general groundsi. 

We have analyzed the simplest gapped excitations. In the repulsive case, when the number of electrons is held 
fixed, the spectrum consists of two gapless holons, as well as a gapful singlet residing in the upper Hubbard band. We 
have calculated its dispersion for various values of the filling and interaction strength, see figure (1,2,3). The gapped 
excitations are independent of the concomitant holon excitations away from half-filling. A dual picture holds for the 
attractive case. 

The upper Hubbard excitation we discussed was based on a two-particle bound state; other excitations emerge 
when higher composites are studied, and can be identified with elementary m-complexes, see eq.(|37]). These have 
higher gaps with respect to the ground state and constitute the upper bands. 



IV. SUMMARY 



We have studied in this article the gapful excitations of the one-dimensional Hubbard model, constituting the upper 
Hubbard band. We have introduced a simple and intuitive construction to incorporate complex momenta solutions 
that correspond to these states. Our analysis has the following results: 

1) In infinite volume states with complex momenta are (elementary) bound states of the hamiltonian, formed of 
any (even) number of electrons. These states are spin singlets and live in the upper Hubbard band for U > and the 
lower band for U < 0. They are renormalized by interacting with other particles, bound or unbound. 

2) These states exist for finite volume if Composite Boundary Conditions (CBC) are introduced. Their existence 
can be proven exactly without postulating a string- hypothesis, via the Bethe Ansatz for Composites (BAC) approach. 
We have shown that all states can be obtained within this scheme. 

3) We clarified the physical interpretation of the SU (2) c ^ arge symmetry in the framework of BAC and have shown 
that it corresponds to the addition of completely local bound states. A simple proof for the corresponding lowest 
weight property of the Bethe states follows. 

4) We have argued that the states with complex momenta cannot have the string form for large but finite volume 
and periodic boundary conditions (PBC). That means that the k — A string hypothesis may give correct results in 
the thermodynamic limit but fails for finite volume and PBC. 

5) The simplest gapped elementary excitations involve a multi-particle bound state. When the number of electrons 
is held fixed the presence of such a state is accompanied by the appearance of at least two holons in the repulsive case 
or two dressed electrons in the attractive case. The gapped excitation is an independent mode only away from half 
filling, becoming a non-dispersive gap at half filling. 

We conclude with two conjectures. The first is related to the possibility of the BAC construction, which follows 
from the composite boundary conditions: we conjecture that the CBC incorporate a symmetry algebra which should 
be similar to the Yangian symmetry found up to now only in infinite systems. 

Furthermore, we conjecture that for a finite system with periodic boundary conditions the Hubbard model is no 
longer integrable in the strong sense; that means the wavefunctions in the iV-particle sector with spin S cannot be 
parameterized by the set {kj, A 7 } of N momenta and N/2 — S spin rapidities alone, if some of the momenta are 
complex. 

We wish to thank A. Ruckenstein, T. Kopp, R. Fresard, C. Bolech, A. Jerez and P. Zinn - Justin for carefully 
reading the manuscript and for their insightful comments. 



A. Appendix - Derivation of eq. (pa) 



This appendix contains the detailed derivation of (|25|), thereby showing the stability of the two-particle bound state 

ib — 

(23) 



if it satisfies the pole condition (|T^) . We define S^L^ by, 



A [231] ~ D 1(23)V K ' 1 '9 ^[123] V iJ -'J 
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i.e. it takes particle 1 with the real momentum k from the left side of the bound pair to the right side, without 
changing their momenta. In both the initial and final amplitudes, we have particle 2 on the left of particle 3, which 

ub 

1(23) 



means k-i—q and k$ = q + . Therefore, S^ 2 S) is given as the product, 



S$3)(*, = Sis(k, q + )S^(K <T) (118) 

which is equation (p4|). To proceed we write S uu (k a , kb) in the form 

S2(k a , k b ) = \{l + s ab )l + 1(1 - s ab )P ab (119) 

where P ab is the spin exchange operator between particles a and b. The phase s ab depends on k a and k b as, 

sin k a — sin k b — 1% , 

Sab = . r — .1 - (120 

sm fe a — sm kb + i if 



We write 



Then 



sin fc — sin q_ — i% sin fc — sin q + — i~ 

_ 2. s+ — H 2. (121) 

sin fc — sin g_ + iU sin fc — sin q + + i% 



1(23) - I \^ + *+) + \( l - S +)P^ Q(l + 8-) + S-)P 

= \ [(1 + »+)(l + «-) + (1 + s-)(l + s+)P 13 (122) 
+ (1 + s+)(l + s_)P 12 + (1 - s+)(l - s-)Pi 3 Pi 2 ] ■ 

Now we use the fact that the spin space of the three particles is restricted: 2 and 3 are in a mutual singlet in 
region [123]. The spin state space in this region is therefore Va = V\ <8 y^ nglet and a two-dimensional subspace of 
the eight-dimensional spin space of the three particles. Under this condition -P 2 3^4[i23] = — ^-[1231- It follows then 
-Pi 3 -Pi 2 = P12-P23 = — Pi 2 if acting on Va- A further identity valid if the operators act on Va reads P12 = 1 P13. 
Note that these operators do not leave Va invariant. In general, therefore, the singlet state of particles 2 and 3 will 
be destroyed upon scattering with particle 1. It is due to a non trivial cancelation of terms if momenta q + , q~ satisfy 
the pole condition for S23, that Sffl 2 3) indeed leaves Va invariant and acts as a pure spin independent phase on the 
wavefunction. We use the identities above to simplify expression ( |122| ) and find 

Si U (23) = ~[(1 + «-)(! + *+) + 2*+(l - «-)] + Pi 3 [l + «- - 3s + + 8+sJ\. (123) 



Explicit calculation shows that the term multiplying P13 vanishes for q + , q satisfying the pole condition (|l^) . It 

2ub 
'l(23) 



follows that S^ns is indeed proportional to the identity on Va, leaving the bound state invariant up to a phase: 



™6 1/-, , . ^ sin fc- 0(g) 



^3)=^ + g -) = sinfc _^ + .| (124) 

In an analogous manner one writes for the S-matrix of two bound states, consisting of particles 1 and 2 parameterized 
by 4>(qi2) and 3 and 4, parameterized by 0(934), 



and one finds immediately 



ebb aub nub /-i nc\ 

D (12)(34) — J 1(34) D 2(34) l lz0 J 



^(12)(34) - T~ I — • U^o; 



f>34 + i 
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= 1. 



The latter is derived as follows. Take the product of equations 



over all j with the product of equations ( |6lj ) over all I: 



nnn 



A-, 



sin kj — if A; — sin kj 



if A 7 - A, 



A- 



sin kj + ij Ai 



sin kj + if Xj 



Ai+ii 



' 1 2 c «£j 



Using now eq.(^) and the trivial identity, Y\ M ^ 



. — = 1, it follows, 

N b AI" 



A 7 — A; — 1% f-r Ai — sm fc, 



l 7 



A-, - A; + if 11 A; - sin kj + if 



(127) 



Equations (|6o|), which set the right hand side of eqn (127) equal to 1, have as a consequence, as asserted, that the total 
momentum is quantized. ■— ■ 
It was attempted by Essler et alEj to show that eq.(|60|) follows algebraically from the set ( p r 8|j5s| , |6l| ) . They restricted their 

attention |te the simpler case of a single bound state and set e^- 4 ! k j+ 2 i) L _ wa y they simply assumed what was to 

be provencj. 
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